Interaction between superconducting vortices and a Bloch wall in ferrite garnet films 
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A theoretical model for how Bloch walls occurring in in-plane magnetized ferrite garnet films can 
serve as efficient magnetic micro-manipulators is presented. As example, the walls' interaction with 
Abrikosov vortices in a superconductor in close contact with a garnet film is analyzed within the 
London approximation. The model explains how vortices are attracted to such walls, and excellent 
quantitative agreement is obtained for the resulting peaked flux profile determined experimentally 
in NbSe2 using high-resolution magneto-optical imaging of vortices. In particular, this model, when 
generalized to include charged magnetic walls, explains the counter-intuitive attraction observed 
between vortices and a Bloch wall of opposite polarity. 
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Microscopic magnetic potentials offer an efficient and 
often indispensable way to manipulate various tiny ob- 
jects, e.g., to trap and guide Bose- Einstein condensates 
and degenerate Fermi gases of ultracold atoms^^ or to 
stretch and twist DNA strands attached to paramagnetic 
beads^. How precisely one can position the trapped ob- 
ject and how large forces can be applied are determined 
by the steepness of the magnetic potential produced by 
the manipulating device. Typically, an assembly of mi- 
crofabricated wires or coils can generate fields with gra- 
dients up to 10 2 -10 3 T/mi2iSi2. Recently, magnetic ma- 
nipulators were created using fcrrimagnetic films, where 
domain walls create locally even stronger field gradients. 
In particular, ferrite-garnet films (FGFs) were used to 
trap ultracold neutral atoms&, assemble and guide col- 
loidal particles&i^, and manipulate vortices in supercon- 
ductors^. Using in-plane magnetized FGFs, where do- 
mains are separated by Bloch walls, has two strong ad- 
vantages: (i) it is easy to move the domain wall and 
thereby tune the magnetic potential, and (ii) one can si- 
multaneously observe the motion of the manipulated ob- 
jects. This "see what you do" ability stems from the large 
Faraday effect in the FGFs, which today are extensively 
used as magneto-optical imaging (MOI) sensors^. In an 
optical polarizing microscope configuration the FGF al- 
lows direct visualization of both the stray field from the 
manipulated magnetic objects and the Faraday rotation 
in the wall itself. 

In this work we present a theoretical model for how 
Bloch walls can function as magnetic manipulators, using 
as example their interaction with a lattice of Abrikosov 
vortices in a type-II superconductor. The configuration 
considered is that of a FGF located near, but a finite 
distance from the surface of a flat superconductor. It 
is shown that the model, generalized to include charged 
walls, explains how the vortices are attracted to such 
walls, and predicts quantitatively the non-uniform flux 
density distribution we find experimentally using MOI. 

Fig. [U shows an image of the magnetic field distribu- 
tion near two linear Bloch wall segments, which them- 
selves appear dark in the image. The superconductor 
is a sine le crystal of NbSe^ cooled to T = 4 K m a 
0.3 mT perpendicular magnetic field, which transformed 




FIG. f : Magneto-optical image of the vortex distribution near 
a Bloch wall in a Bi-substituted lutetium iron garnet film 
placed on top of a superconducting NbSe2 crystal with tran- 
sition temperature of 7.2 K. The slightly uncrossed analyzer 
and polarizer setting used here implies that the dark wall and 
the bright vortices have opposite field polarities. Image di- 
mensions are 70x70 /jm 2 . 



into quantized vortices when entering the superconduct- 
ing state. Each vortex is here seen as a bright dot. Ev- 
idently, the magnetic wall has a considerable attraction 
on the vortices, since their number density increases near 
the wall. This observation forms the experimental basis 
for the theoretical modeling presented below. Interest- 
ingly, a most surprising feature visible from Fig. [T] makes 
the problem even more challenging. There is opposite 
contrast between the dark wall and the bright vortices, 
which means they are of opposite magnetic polarity. In 
this case one would expect from the models presented 
previously in the literatur o 14 i 15 ' 16 i 17 that a Bloch wall 
should repel the vortices. As will be shown, by account- 
ing for additional magnetic charges due to misalignment 
between the wall direction and the in-plane magnetiza- 
tion vector within the domains, the sign of the interaction 
can become inverted. 

The two Bloch wall segments seen in Fig. [1] are actu- 
ally part of a larger zigzag pattern. Extended zigzag do- 
main walls are commonly present in FGFs with strong in- 
plane anisotropy-i£. An example is shown in Fig. [2] (top), 
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FIG. 2: Top: MO image showing a zigzag Bloch wall in a 
FGF separating two domains with antiparallel in-plane mag- 
netization. Bottom: Sketch of an in-plane magnetized FGF 
with a Bloch wall placed above a superconductor. The mag- 
netic charges along the vertical sides of the wall can lead to a 
net attraction between a wall and vortices, as seen in Fig. [1] 



where the zigzag line separates two domains with an- 
tiparallel magnetizations that meet head-on. By folding 
into a zigzag pattern the domain boundary reduces the 
density of magnetic charges at the wall, which helps to 
minimize the energy^. To describe the interaction be- 
tween one segment of such a zigzag wall and a super- 
conducting vortex, we introduce the model illustrated in 
Fig. [5] (bottom). The superconductor occupies the half- 
space, z > 0, and the wall is directed along the y axis, 
which forms an angle ip with the magnetization direction 
of the two domains. Experimentally we detect only tiny 
changes of the wall width, 2W, during cooling through 
T c , in agreement with [Til. Thus, in calculations the wall 
is approximated as a fixed, uniform, out-of-plane magne- 
tization M z (|x| < W) = —M s . Inside the domains there 
is an in-plane magnetization with a component normal 
to the wall given by M x (|a;| > W) = —M s (x/\x\) smip, 
where M s is the saturation magnetization. For ip = 0° 
the present model reduces to the non-charged wall case. 
The My component is omitted in the analysis since the 
wall is assumed to be infinitely long. 

Stray fields from the wall induce shielding currents in 
the superconductor, which we determine using the Lon- 
don theory. The equations valid inside and outside the 
superconductor then read 



\- 2 A-V 2 A =0 



-V 2 A = ^o(V x M)j 
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z < 
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where A is the London penetration depth and the vector 
potential is A = Ay. The shielding currents flow in the y 
direction and their density equals J y — —A/[Iq\ 2 . A vor- 
tex present in the superconductor then feels two forces. 
First, the direct force from the FGF, which can be cal- 




FIG. 3: The calculated forces on the vortex from the Bloch 
wall: F ± is repulsive and F" is attractive. Their sum F vw 
changes sign at x* ~ 1 /im. The parameters used here are: 
sin if = 0.34, 2W = 0.6 /im, h = 0.8 /im, a = 140 nm, 
A = 200 nm, and M s = 50 kA/m. 



culated from the free energy term, /iq J M ■ H v dV , where 
H v is the stray field from the vortex. Second, the Lorentz 
force from the shielding currents in the superconductor, 
= J(r) x <J> integrated over the length of the vortex. 
$0 is the magnetic flux quantum, and we will simplify 
the treatment by assuming the vortex to be straight and 
aligned with the z axis. Interestingly, the two forces turn 
out to have exactly the same magnitude and direction, as 
was noted also in Ref. [HI where a similar configuration 
was analyzed. 

It is convenient to express the total force on the vortex 
in the x direction as 



pvw _ jpj 



F" 
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where F 1 - and F" are the contributions from the per- 
pendicular and in-plane components of M, respectively. 
Their Fourier transforms arc obtained as 
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where r = \J A -2 + fc 2 , a is the gap between supercon- 
ductor and FGF, and h is the FGF thickness. For the 
configuration illustrated in Fig. [5] the force is always 
attractive, whereas F 1 - is repulsive. This qualitative re- 
sult can be easily understood by considering the interac- 
tion between the magnetic charges involved. The stray 
field from a vortex is closely approximated by that of a 
magnetic monopole located at z <~ A and with strength 
— 2$cr^- Thus, the vortex is attracted to the positive 
charges along the vertical sides of the wall and repelled 
by the perpendicular dipole charges. The charge repre- 
sentation also yields the correct magnitude of forces given 
by Eqs. ([3]) and ((4]) in the limit A — > 0. The supercon- 
ductor then perfectly screens the magnetic field created 
by M and its presence should be accounted for by adding 
the corresponding mirror charges. 
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FIG. 4: The excess vortex density 8N(x) near the Bloch wall 
for various wall width 2W calculated using L — 160 /im 
and the other parameter values as indicated in the caption 
of Fig. E 

From the inverse transform of and we ob- 
tain the spatial variation of the two force contributions, 
which arc plotted in Fig. [3] together with the total force 
on the vortex. At sufficiently large x the magnitude of 
F'I is larger than that of F ± , which is expected since 
the monopole-monopole interaction should dominate at 
long distances. However, at short distances F 1 - becomes 
dominant, and the total force changes sign at some dis- 
tance x* . For our set of parameters, the repulsive region 
| a; | < x* is very small, less than 1 /im. This is why the 
repulsive region under the wall is not visible in the image 
of Fig. [TJ This also explains why we observe the counter- 
intuitive attraction between the Bloch wall and the vor- 
tices of opposite polarity. A related phenomenon when 
ferromagnetic domain wall stimulates superconductivity 
due to its stray fields was considered in Ref. [2ll . 

We consider next an initial state with a uniform dis- 
tribution of vortices in the superconductor, and a sub- 
sequent introduction of a Bloch wall. This results in a 
perturbation of the vortex density, 5N(x), which creates 
an additional force acting on every vortex. In the equilib- 
rium, the additional force everywhere balances the force 
from the Bloch wall, i.e., 

where F vv is the x component of the repulsive vortex- 
vortex interaction, and the vortex matter is considered 
as a continuum. This equation represents a perfect 
shielding of the domain wall by the vortex matter. In 
Fourier space, the perturbed vortex density becomes 
6Nk = F£ w /F£ v . The vortex- vortex interaction can 
be obtained from the currents around a flux line in a 
half-space, first calculated by Pearl^ 2 ., and gives 

J™ = ^ \ ik (\k\L + ~—^—) (6) 
k A 2 \k\r 2 V A 2 t |fc| +tJ 

where L is the flux line length. The term proportional 
to L is the conventional (Abrikosov) bulk contribution, 
while the other term is the surface contribution. 
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FIG. 5: Top: Distribution of vortices formed in the pres- 
ence of a Bloch wall. The image was taken after the wall, 
seen in Fig. [T] was removed by an in-plane field of the order 
of a few /J.T applied perpendicular to the indicated x-axis. 
Bottom: Vortex density obtained from the image (each sym- 
bol represents one vortex) together with the theoretical curve 
calculated for L — 160 [im, $oiVo = 0.3 mT and remaining 
parameters as listed in the caption of Fig. [3] 

The resulting perturbation 8N(x) of the vortex density 
induced by a Bloch wall is shown in Fig. 2J Its profile is 
strongly dependent on the wall width 2W. In the small 
W limit the total force on a vortex is everywhere attrac- 
tive, hence the vortex density increases monotonously as 
one approaches the wall. For increasing W the density 
profile develops a minimum below the center of the wall, 
and the resulting 6N(x) becomes non-monotonous with 
a minimum at the center and maxima near the two wall 
edges. For sufficiently large W the vortices right below 
the wall are expelled creating a narrow depleted area. 

The theoretical vortex density profile 8N{x) can be 
compared to our MOI observations of the vortex distri- 
bution near the Bloch wall. The FGF had a thickness h = 
0.8 /xm, and saturation magnetization M s = 50 kA/m. 
The film, with no additional layers (contrary to stan- 
dard MOI indicators), was placed directly on top of the 
0.3 mm thick NbSe2 crystal with a gap of a = 140 nm, as 
determined from the optical interference patter 1*22. This 
gap equals one quarter of a wavelength, and gives optimal 
transmission. 

The vortices were formed slightly below T c where flux 
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pinning is negligible. Thus, one expects that the vor- 
tex positions seen in Fig. [1] represent a frozen picture of 
an arrangement where the vortices adjust only to bal- 
ance the interaction with the wall. To obtain a better 
view, we made use of the fact that further cooling to 4 K 
increased the vortex pinning considerably, and removed 
the Bloch wall without creating noticeable change in the 
vortex positions, see and Fig. (top). From this image, 
the positions of all the vortices inside the marked rect- 
angular area were identified, and the Wigner-Seitz cell 
of each vortex was determined using standard triangu- 
lation^i. The local vortex density was obtained by in- 
verting the cell area, and is shown in Fig. [5] (bottom) for 
every vortex versus its coordinate x. The vortex density 
near the wall clearly exceeds the background density of 
$o^o ~ 0.3 mT. The theoretical curve N + 6N(x), plot- 
ted in Fig. [5] (bottom), was calculated using (p = 20° de- 
termined from Fig.[TJ and a wall width of 2W = 0.6 fim. 
The calculated curve reproduces very well not only the 
sign, but also the magnitude of the excess vortex density. 
This agreement was achieved using the penetration depth 
A = 200 run as an adjustable parameter. This value of 
A corresponds to the temperature of 6.8 K (slightly be- 
low T c = 7.2 K) which is thus the temperature when the 
vortices got frozen. The vortex length L=160 /im was 
another free parameter. It is smaller than the crystal 
thickness 300 /im to compensate for the overestimation 
of the Abrikosov interaction term in the continuum ap- 
proximation. 

An open question remains regarding the large appar- 
ent width of the Bloch wall, approximately 3 /im as 
seen from the image in Fig. [TJ The theoretical es- 



timate obtained by minimizing the sum of exchange, 
anisotropy and magnetostatic energies is given implicitly 
by the equation 1 ^ (1 + w)~ 2 + aw~ 2 = 1 - K u //j, M 2 , 
where w = 2W/h is the normalized wall width and 
a = 2TT 2 A/iioM 2 h 2 . Substituting the effective exchange 
constant A = 2 x 10~ n J/m and the uniaxial anisotropy 
constant K u ~ 10 3 J/m 3 [l| we obtain 2W w 0.6 fxm. 
The discrepancy between the observed and estimated 
wall width is probably due to the optical diffraction which 
significantly distorts the image of objects of the order of 
the light wavelength, 0.55 /Ltm. 

In conclusion, mobile domain walls found in in-plane 
magnetized ferrite garnet films were investigated for pos- 
sible use as magnetic micro-manipulators. It was shown, 
choosing superconducting vortices as a case example, 
that such films can serve to both apply forces and si- 
multaneously monitor the results of the action. A theo- 
retical model for the interaction was developed, with the 
vortices described within the London approximation, and 
the domain wall represented by a charged magnetic wall. 
The charged wall model, which includes magnetic poles 
on all the sides of the wall's rectangular cross-section, 
is shown to give a very good quantitative description of 
the attraction of vortices to such a wall. The compari- 
son was made by direct observation of individual vortices 
using the magneto-optical imaging technique. 
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